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RE/XTIONS  OF  SII.3PIE  SYSTEIS  UiffiER  BL/iST  LOADING 


by  D.  Ilontgomery  and  A.  H.  Taub 


Abstract 

The  differential  equation  I'.6c+ F(x)  =  p(t)  is  considered  for  seme 
simple  cases  of  blast  loading.  The  right-hand  side  is  assumed  linear, 
and  F(x)  on  the  one  hand  is  taken  as  constant  and  on  the  other  is 
taken  as  linear  fr(»n  the  origin  to  the  constant  and  then  as  remaining 
constant  for  larger  values  of  x.  It  is  shown  that  the  situations  in 
the  tvro  cases  differ  moderately.  An  approximation  formula  is  de¬ 
veloped  by  Ti-hich  certain  information  in  the  latter  case  can  be  ob¬ 
tained  from  the  former. 

1 .  Introduction 

In  discussing  the  behavior  of  various  targets  under  blast  loading  it  is 
often  possible  to  reduce  the  mathematical  problem  to  that  of  a  one -dimension¬ 
al  system  governed  by  the  equation 

M  ^  +  F(x)  =  p(t),  (1) 

dt*- 

Trfiere  x  is  the  displacement  of  thq  system,  F(x)  is  the  restoring  force,  M  is 
the  equivalent  mass  of  tlic  system,  and  p(t)  is  the  force  [  *  pressvire  x  area] 
acting  on  the  system  where  .p(t)  is  dependent  on  time. 

Equations  of  this  form  arise  in  many  problems?  for  example,  if  the  tar¬ 
get  is  elastic  and  has  various  modes  of  vibration,  its  response  is  determined 
by  solving  a  set  of  equations  .of  the  type  of  Eq.  (l )  v;hcre  F(x)  is  of  the 
form  Again,  this  equation  is  found  in  the  treatment  given  by  Christo- 

pherson-5-^  in  R.C.  3h9  of  the  action  of  brick  T.alls.  There  it  is  shovm  that 
P(x)  may  be  replaced  by  a  constant. 

In  the  application  we  hsave  in  mind  (blast  vjave)  the  function  p(t)  is 
zero  for  negative  time,  has  a  finite  initial  value  p^  at  t = 0,  decreases  to 

1/  "A  modification  of  the 'impulse  criterion  for  blast  damage,”  by 
D.  G.  Christopher son,  R.C.  3h9f  Sept,  ^9h2  (Confidential), 
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zero  again  at  time  tQ,  and  becomes  negative  thereafter*,  rising  to  zero  at 
sQae  •  later  time .  The  problem  Tri.th  rrhidh  ire  are  mainly  concerned  nay  be 
stated  as  follows: 

Tiihat  relation  must  exist  between  and  t  in  order  that  the  maxi¬ 
mum  of  the  solution  of  Eq,  (l)  be  a’  speciiied  quantity,  say  X3? 

If  tlic  solution  of  this  problem  is  kncavn,  then  for  a  target  such  as  a  brick 
wall  we  can  determine  the  relation  between  p^  and  t^  that  Trill  just  cause 
the  target  to  fail  —  that  is,  reach  a  critical  displacement  with  zero  veloc¬ 
ity.  The  quantity  p^  is  called  the  peak  pressure  acting  cn  the  target,  and 
t^  is  called  the  dtiration  of  the  pressure  wave.  The  area  under  the  pressure¬ 
time  curve  betrreen  t  =  0  and  t  =  t^,  called  the  positive  impulse,  may  be  re¬ 
lated  to  Pq  and  t^.  The  result  may  then  be  e:q>ressGd  in  terms  of  the  peak 
pressure  and  pqsitive  impulse  just  necessary  to  cause  failure.  If  the  rela¬ 
tion  betTjGon  peak  pressuiro  and  impulso_acting  on  the  target  and  the  same 
quantities  in  the  blast  mve  are  knoim,  then  for  any  charge  vfcight  a  dis¬ 
tance  can  be  determined  that  is  the  limiting  distance  at  which  the  target  is 
destroyed.  In  order  to  perform  the  last  calculation  the  dopondonce  of  peak 
pressure  and  impulse  in  the  blast  wave  on  woi^t  of  charge  and  distance  must 
boJcno\'nl.  These  quantities  miost  bo  corrected’  for.  reflection,  diffraction, 
and  motion  of  the -target  in 'order  to  o’btain.the  peak  pressure  and  impulse, 
acting  on  the  target.  -  ■  '  ■  -• 


This  paper  will  be  concerned  ivith  the  determination  of  the  relation  be 
tvjDon  peak  pressure  and  impulse  acting  on'  the  target  for  a  given' maximum 


displacement  for  special  cases  of  Eq 
f  olloiTS : 

A:  p(t) 

Case  I:  F(x) 
or 

Case  II:  F(x) 


(l).  The  specializations  made  arc  as 


=  constant  =  P, 


.'■p  ■  ■  ■  . .  •  • 

:;q- X,  'O^x  ^XiS  (a) 

■^1 

< 

P,  X  >  Xi.  (b) 
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Case  I  is  a  limiting  case  of  case  II.,  If  the  desired  deflection  is  X3 
and  if  x^/xs  approaches  zero,  then  general  existence  theorems ,  gtiarantoe  that 
the  solutions  in  case  II  approach  those  in  case  I.  Hovrever,  there  is  no 
guarantee  that  a  given  value  of  x^/xs,  say  0,01  for  example,  vd-ll  bring  the 
solutions  near  each  other  to  an  accxiracy  of  about  the  sane  size.  Actually 
•JO  find  that  the  solutions  can  differ  to  a  greater  degree  than  0.01  in  this 
'  case,  although  the  difference  is  not  excessive.  Wo  exhibit  numerical  calcu¬ 
lations  bearing  on  this  point,  and  vra  also  develop  a  formula  that  makes  it 
easy  to  calculate  from  the  limiting  case  tdiat  the  situation  is  for  a  given 
value  of  X1/X3  provided  this  value  is  not  too  large.  We  consider  only  cases 
T;here  X3>  x^  since  in  such  cases  a  target  vdll  be  destroyed  idien  it  reaches 
a  deflection  X3  'Tith  zero  velocity.  At  the  end  i7e  also  take  up  a  related 
question  vdiosc  description  wo  postpone. 


2.  Solution  for  case  II 


Wo  shall  treat  case  I  as  a  special  case  of  case  II  and  proceed  first  to 
obtain  the  solutions  in  the  latter  case. 


In  the  interval  from  0  to  x^  the  solution  is  as  fdllori’ra: 
..  .  Po 


M<.? 


sinut  t  j.  ^  ,  j.  j 

^  V"  ^  ' 


■vdicrc  <j^  =  P/ifei,  and  hence  in  this  interval 


X  = 


Po 


[cos  (Jt  +  wt-  sin  wt  -  1  ] . 


llw^t  "  "o 


(2) 


(3) 


Let  ti  be  the  timc_at  v.Mch  the  displacement  reaches  x^.  liiking  use  of 
the  fact  that  =  P/litc^,  -kTc  see  that 


X  -  ^ 
p 


sxn  £ot  X  t  T  ,  ^  , 

- 7 - ^  +  1  -  cos  <jti 


-  °  -1  . 

Dividing  by  Xj.  and  rearranging,  vro  find  as  the  equation  determining  t^ 

P 


wt^  cos  cjt^+  wt^  -  sinot^ 


ut„  1  - 


ih) 
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Tie  let  Xj  be  the  j  value  of  x  at  ti  and  'ive  denote  by  I  the  quantity 
•Ip^t^ '.Thich  is  the  area  under  the  curve  p(t)  from  0  to  t^.  Then  from  Eq: 
(3) 

*T^  “  tj2^  sin  wti  +  cos  tjti  -  1). 
o 

Tfiien  t  is  greater  than  ti  and  x  is.  greater  than  Eq.  (l)  be  cones 


(5) 


=  p(t)  -  P’  =  (Po  *■  ^ 


I^kking  use  of  the  fact  that  x  =  Xi  and  x  =Xi  nhen  t  =  ti,  ttd  find  that 

'  (■]  .-  ii'i  »  -D  ('i  j.  ^  Po  fj.  j. 


lik  =  i.&^  +  i^Poh  (t  - 1^)  -  (t  - 1^) 


(6) 


and 


f  * 

Ibr  =  I&^+I^^(t-t^)  +  J-^p^(l-|ij-Pj  (t-.t^)^-|^Jt-t^)f  (7) 

Ijet  ta  be  the  value  of  t  at  t^iich  the  solution  given'  by  Eq.  (7)  has  its 
maximum  value,  and  let  X3  be  this  maximum.  T/hen  t^ta  the  left-hand  side  of 
Eq.  (6)  is  zero  and  vre  obtain  '  ■ 

^  “  ,po  (‘  "*  ^)'  • 

bet  T3  =  (ta- ti)/tQ.  Rearrangement  gives  r. 


1.6c 


1  =  mm  2 


■(r-ii)-il 

\  ^ol  Fo 


Solving  for 


T3  =  (1  -  b;  _  +  (1  _  -  -I:\ 

"  V  Poj-  I  .  I'  Vo) 


.1/2- 


(8) 


idierc  Tve  must  choose  the  "si^  before  the  ^qviare  root  so  as  to  make  fa  posi¬ 
tive.  In  developing  the  apprcodlmation  formula  itg  consider  the  case  iThere 
1  -  tyto-  P/Po  is  positive. 

To  find  Xa  'CJe  substitute  this  value  in  Eq.  (7); 


I6C3  "■  Ii6c3_  +  toT  3 


I.XK1  +  I 


hk)- 


+  r  -  rz 

^o  *3  g  3  r. 
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Replace  by  its  value  from  the  equation  preceding  Eq.  (8) : 


I&s  *  I.bc 
=  £2  p  =  ££  1, 


.  Po^o^a 
1  ^  — 


2T3 -  3  1  - 


Novr  Pq  ~  ^  ^  ~  ^  and  hence 


-0 


(8a) 


(9) 


Instead^of  plotting  p^  against  t^  it  is  more  convenient  to  plot  l/VSPMxa 
against  Pq/Fj  a^nd  vie  shall  next  derive  a  fontiula  expressing  the  first  of 
these  quantities  in  terms  of  the  second  in  the  limiting  case.  In  this  case 
ti  =  Xi=0.  Here  Eq.  (9)  becomes  meaningless  because  in  Eq.  (9)  vre  have  -used 
the  expression  P/xi  for  a  slope.  Kovrever,  fran  Eq,  (8a),  in  this  case 

>*..  -  ^o^IPo 

11X3  -  - T - 


2r3-  3(1 

'-F 

L  \ 

Pq/  ^ 

and  also  in  this  sane  case  by  Eq.  (8) 


r,  =2(1 


Hence 


and 


2PltC3  = 

3Pn 


I2  ^  3  Po 
2Pi.&  3  TS  P 


-0 


FF 


(10) 


(11) 


From  Eq.  (11 )  a  table  of  values  may 
be  computed  relating  PQ/P  and  l/v/2PItc3 . 
Table  I  and  the  graph  of  X3yx3  =  0  in 
Fig.  1  present  these  values. 

Notice  that  for  some  computations  in 
case  II  it  is  convenient  to  use  the  fol- 
lo\-jing  relation 

1  £0 

I  _  2 


Table  I .  Values  of  l//2PI.bC3  in 
the  limiting  case . 


P  _  1  Po 


V'2I.'Fu)2xiX3  P  »X3 


(12) 


Po 

I 

Po 

'  I 

,p 

\^PMX3 

‘F 

v^Plbca 

1.2^ 

5.i;13 

6. 

1.391; 

1.^ 

2.756 

7 

l.UIii; 

2 

1.732 

8 

I.U96 

1.378 

9 

1.550 

h 

1.333 

10 

1 .601; 

5 

1.353 
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The  quantity  I^/2P3fe3  has  the  foUoiving  physical  interpretation;  it  is  the 
ratio’  of  the  kinetic  energy  given  to  the  target  if  the  loading  is  truly  im¬ 
pulsive  (the  impixlse  I  is  communicated  before  any  displacement  or  velocity  is 
acquired  by  the  target)  to  the  static  -OTrl^:  done  pn  the  target  rriien  it  is  dis¬ 
placed  to  failure.  This  ratio  v;ould  bg  one  for  impiolsive  loading.  Actually 
in  the  limiting  case  this  ratio  is  a  function  of  P/Pq  and  its  minimum  value 
is  U/3.  Thus  the  fact  that  the  loading  is  spread  out  over  a  finite  time  has 
an  appreciable  effect  on  the 'behavior  of  the  system. 

The  fact  that  the  value  of  rises  slowly  for  values  of  p^/P 

greater  than  four,  junplies  that_ih  this  range  of  P^/P  the  ’’impulse  criterion'! 

I 

is  approximately  true.  That  is,  if  the  impulse  in  the  pi’essure  rave  acting  <x 
the  structure  is  greater  than  and  approximately  the  minimum  value  the  target 
Tiill  break,  provided  of  course  p^/P  is  greater  tlian  foiu:.  In  the  range  Tdiere 
Pq/P  is  less  than  tiTO,  but  greater-  than  one,  the  value  of  l/v&Pilxs  changes 
very  rapidly  for  si3.ll  changes  in  p^/P.  This  means  that  the  brealdng  of  the 
target  is  folloi'dng  a  pressure  criterion. 

3*  Approximation  formulas 

Tfc  shall  noiT  find  a  method  to  obtain  approjcimately  the  value  of  l/vSPLtxs 
for  a  given  value  of  Xx/xa  from  the  value  of  l/v^PI-ks  in_tho  limiting  case. 
Tlio  eqijations  we  need  for  this  purpose  are  Eqs.  (U),  (5),  (8),’ and  (9). 

For  convenience  we  collect  these  foimiula-s  in  one  place 


tjt^  cos  ut  +  ot  -  sin  ut  =  «t„(l - 

O  17.  1  o  (  p 


Itc  2 

-  -  ■  (totg  sin  «Jt-  +  cos  <jt-  -  l), 
t^ 


Wi  .  .  zf 

\  *0  Po/ 

I 

!- 

0 

I _ 

1/2- 


2^3- 3(1  -h-.£ 


Po 


ih) 

(5) 

# 

(8) 

(9) 
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Suppose  that  P/Pq  is  fixed.  If  vre  also  fix  t^t^,  then  tot^  is  deter¬ 
mined  by  3q.  ih)  •  Then  l&i/l  is  given  by  Eq»  (5)  and,  since  ti/to  wti/utQ, 
I'a  is  given  by  Ec.  (8)  and  x^/xs  is . determined  by  Eq.  (9).  Hence  for  each 
value  of  P/Pq  there  '.ull  be  a  value  of  cot^  vihich  gives  X^/xg  a  fixed  value. 

Assuming  that  x^/xg  and  P/pq  are  fixed  we  now  estimate  \’iiat  the  A?alue 
of  t>tg  vrLH  be  -when  x^/x^  is  small. 

For  a  first  rough  estimate  assume  that  ti  =  0  and  that  T3=2(1-P/pq) 
as  in  the  limiting  case.  Then  from  Ec.  (9) 


£3 


“if’--)- 

°\  Po/ 


■1  J  \ 

IVe  may  drop  the  one  as  unimportant  rfaen  Xa/xi  is  large  and  get  as  a  first 
estimate  (otQ)^  of  wtQ, 


Uo)l  = 


X3  P 
^  ^  Po 


2  1 


--•r 

■  Po'/ 


(13) 


Tfe  nov.'  assume  that  a  second  approximation  (oJtQ)2  is  given  by 

'(<jto)2  =  /3(tjto)^, 

where  ^  is  a  quantity  to  be  determined. 

From  Eq.  (h)  it  is  seen  that  a  good  apprco:inis.tion  to  cjti  is 

<ut-i  =  arc  cos  (l - ^  |, 

\  ,Po/ 

and  from  Eq.  (5)  an  approximation  for  l&j/l  is 

I'lXi  _  2  sin  wti 
—  ZIt;  • 

Tfe  also  make  the  follcr^zing  estimates  for  T3  and  T3: 

1 


(iu) 


(15) 


(16) 


_  0/  p\  1 

\  1>0  Po/ 


(17) 


and 


^2 


^  -  if-  ^  ^  ^  -  s)- 
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The  valxB  has  been  chosen  so  that  the  desired  value  of  Xa/xi  is 

given  by 


^3  -  2  Po  (.  \Z  ('  P  f  •  • 


(19) 


but  substitution  of  Eq.  (1 8)  in  Eg.  (9)  and  use  of  the  second  approximation 
(‘^1'o^2  “  1  gives 

Equating  these  two  values,  again  dropping  the  one  as  imimportant  ivhen 
Xs/xi  is  large,  and  dividing,  Tre  find 


Hence 


Yjhere 


ot  = 


1 


7TZ 

?o  . 


1  +  CC  ^ 


3/2  ^ 

Po 


Substituting  the  relation  given  by  Eq.  (I6)  and  remembering  that  iiAo 


sin  4oti 


(utQ)^  * 


(1  -  if  . 

\  Po/ 

and  then  using  the  relations  given  in  Eqs.  (13)^  (ih),  and  (l^). 


OC  =  l  3 


73 


[3in{arc  =o=(l  - i)| -(l  .  {l-i)]| 


\  Po/  X, 


Po/ 

•D 
^  O 


1/2 


or 


06  = 


I—- 

7VX3 

A  ' 
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^2  = 


P  = 


1 


1  +  CL 


1  + 


7VX3 


[E 

1x3 


7  V  X3  r  »2xi/x3 
- : — i  +  .11  +  2 — fii — i  , 

2  V  U 


(20) 


Table  II.  .  Values  of  Z  and  fl. 


Po 

.p 

7 

7 

? 

7" 

T 

^  for 

—  =  0.01 
X3 

1.25 

'll. 323 

2.162 

2:337 

0.760  ' 

1.5 

2.767 

I.38L 

0.958 

.'852 

2 

1.677 

0.838 

.351 

■  .920 

3 

0.959 

.L80 

.115 

.953 

h 

.688 

.3hh 

..0592 

.966 

5 

.522 

.261 

.03L1 

.97li 

6 

.L25 

.212 

.0225 

.979 

7 

.359 

.180 

.0162 

.982 

8  • 

.311 

.156 

.0122 

.98L 

9  „ 

.27L 

.137 

.009L 

.986 

10 

.2L5 

.122 

.007L 

1  .988 

Thus  vie  have  achieved  our 
purpose  —  to  find  the  approxi¬ 
mate  value  of  A  table  of 
values  of  7/2  and  7^/8  as  well 
as  the  values  of  ^  when  X1/X3  =0.01 
is  given  in  Table  II. 

If  we  replace  in  Eq.  (12) 
by  (cjtQ)j^  we  see  tl:iat  we  obtain 
the  value  of  l/\/2Pl.fcc3  for  the 
limiting  case.  Hence  the  factor 
^  is  also  the  factor  which  -vhen 
multiplied  by  l/V^PTS^  in  the 
limiting  case  yields  the  value 
(approximately)  for  any  \’alue  of 

X1/X3. 


Table  III  shows  l/vSHbEs  as  computed  in  certain  cases  and  as  given  by 
the  approximation  formula  derived  above.  It  can  be  seen  that  the  approxima¬ 
tion  formula  is  quite  accurate  in  these  cases. 
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Hence,  to  get  a  good  esti¬ 
mate  of  l/v^P%3  for  a  given  small 
value  of  .Xi/xs,  ccanpute  ^  from 
Sq.  (20),  using  ip  many  cases  the 
vallies  of  V/2  and  ?'®/6  from  Table 
IT .  Then  multiply  ^  and  the 
.  limiting  value  of  l/v2PI5c^  from 
Table  I.  This  approximation 
formula  is  quite  accurate  vdien 
1/^3  is  small,  and  is  fairly 
accurate  for  values  of  x^/xs  as 
large  as  0.3  or  O.k-  Figure  1 
gives  the  values  of  l/y/ZPlhz^ 
for  various  values  of  x^^/xa  be¬ 
tween  0  and  I  -  .  The  case  Xi/x3  =  1 


Table  III .  Comparisoii  of  exact  and 

'appro:cimate  values  of  l/v^Pik 


Po 

X3 

Exact  Value 

Of  I 

Value  Given  b 
Appr  oximati  on 
Fomula 

p 

1.^ 

0.020^ 

2.230 

2,265 

2.0 

.0h2h 

1.U57 

i.h58 

2.0 

.012li 

1.^77 

1.578 

3.0 

.01^0 

1 .300 

1.299 

3.0 

.0069 

1.32h 

1.32h 

h.O 

.0223 

1 .269 

1 .266 

h.O 

.0129. 

1 .263 

1 .281 

7.0 

.0091 

1  .h20 

1.U20 

10.0 

.0220 

1^575 

1.575 

10.0 

1 

.0125 

1 

1.^83 

1.532 

may  be  handled  as  follo^re.  In  this  case  F(x)  is  a  straij^it  line  and— 

X  =?^Ti 

'  ^  15^^  L  ^■'-o  j ' 

p ' 


/ 


using  the  fact  that  w 


2  = 


T5e  obtain 


arc  tan  <ato 
<jto 


Vife  find  also  in  this  case,  by  Eq.  (12), 


I  =  1  ^o  , 

2V2  p 


The  graphs  for  this  case  and  for  the  cases  utiere  x^/y^s  -  0.0,  0.01,  0.1,  0.2, 
and  0.5  are  shcum  in  Fig.  1 . 

The  curves  of  Fig. .1  all  have  vertical  asymptotes  on  the  left.  To  find 
them  proceed  as  follo^vs.  Loolcing  at  the  Eqs.  (h) ,  (o),  (8),  (9),  and  (12), 


The  follcT.'ing  equation  is  derived  by  solving  tlie  differential  equa¬ 
tion  for  x(t).  Determine  the  time  at  vdiich  the  maximum  is  obtained  from  the 
equation  x(tj.)  =  0.  Substitute  this  value  of  t^  in  the  equation  x(t  1)  ®  X3,.  See 
R.C.  6,  ’*The  design  of  buildings  against  air  attack  (Part  2),"  Ibirch  1939* 
Restricted. 
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let  ijt^  approach  infinity.  Then  ^'3  approaches  zero.  Let  ^  §  1, 

be  fixed,  and  attempt  .to  find  the  corresponding  valu^  of  Xj/xs.  Since 
l-tx/tQ-P/po  is  negative,  the  approxinaticm  for  ^ven  iiy  Eq.  (l?)  is 
to  be  replaced  by  the  expression  cbtained  by  the  choice  of  signs  in  Eq.  (8) 
vihich  Kiakes  ^'3  positive.  This  gives 

^  -1  ^  ' 

^  I  2U  -  t^yt^  -  p/:d'J* 

Since  is  negligible  compared  to  1  -  P/p^,  we  obtain 

_  _  sin  tjti 

^  “  ‘•^o(p/p7- 

Substituting  in  Eq.  (19)  and  replacing  sin^otiby  2P/Pq-  (P/p^)^  vre  obtain 

X3  _  .  1  Po  2P/Po  -  (P/Pq)^  r.v  M 

5^-1  +  H  T  (P/Po"^'T)a"“  ^^3  -  3(1  -  ^/Pq)^- 

Ihiking  use  of  the  fact  that  '?3  is  small  compared  to  3(l-  P/Prt)  I’ra  yarite 


^  =  1  . 


1  2  -  P/Po 

^  TT^p- -  --TF 


3(P/Po  -  1) 


and 


.  1 . 

p/Po  -  1 


^  =  2P/Po  2  ’ 
X3  P/Po  ■ 


This  may  be  -written 

£2  =  1  1  -II 

P  ^  “  2  X3* 

Hence  this  is  the  location  of  the  vertical  asymptotes  for  the  curve  asso¬ 
ciated  with  X1/X3.  VJhen  X3_/x3  =  0  the  asymptote  is  at  1,  and  as  X1/X3  in¬ 
creases  to  1  the  position  of  the  asymptote  shifts  linearly  to  -I-. 

As  y/e  have  seen,  the  corves  all  have  a  vertical  asymptote  given  as 
above.  After  this  they  drop  rather  soon  to  a  rainimm  and  then  rise  gradu¬ 
ally.  The  position  of  the  minimm  .varies  from  about  1.3  for  Xx/x3=  1  to 
li  for  Xx/x3  =  0.  \  ' 


^  ..0 


ll.  Comparison  of  case  II  to  an  elastic  system 

Returning  to  the  differential  equation  Eq.  (l),  itb  disci;iss  next  a  prob¬ 
lem  iNhich  arises  in  connection  vdth  F(x)  as  given  in  case  II  and  as  gL'V’en  in 
still  another  case  called  case  III. 

Tile  function  F(x)  =  kx  ifliere  k=  P/xi  is  the  value  given  in  (a)  of  case  II 
Thus  the  curve  of  case  IH  is  merel3''  a  continuation  of  the  straight  line  viith 
•which  the  curve  in  case  II  begins .  Suppose  that  the  desired  maximimi  deflec¬ 
tion  in  case  II  is  X3  and 'that  the  desired  maximum  deflection  in  case  III  is 
x^.  The  area  under- the  curve  II  from  0  to  X3  is 

,  '  loci'  • 

As  “  kXiX3 - 

The  area  under  the  cur've  III  frcm  0  to  X3  is 


A 


3 


Under  some  conditions  it  is  reasonable  to  suppose  that  if  As  “  A3,  then 
the  Pq  and  t^  TJhich  produce  a  maximum  deflection  X3  in  case  II  v.lll  produce  a 
maximum  deflection  X3  in  case  III.  This  conjecture  vdll  be  examined  belorj", 
and  it  •';dll  be  shovm  that  .it  is  not  alvra,ys  true. 


For  Ag  “As,  the  folloiving  must  hold 


x^^  =  2x1X3  -  Xi, 


or  vjhcn  Xi  is  small. 


X3^  “  2x1X3, 


(21) 


Let  Xjxj  be  the  maximum  deflection  in  case  III.  Then 


2p_  ^  arc  tan  ut. 


"ni  , 


1  - 


CJ'bl 


Po 

-p-  Xi 


1  - 


arc  tan  ot. 


U'Cq 


For  small  values  of  Xi  this  is  approximated  fairly  v.’cll  by 


"III 


2 


II  T..n  7 


1 


In  case  II  assume  that  Xi/xs  is  so  small  that  the ’ maximum  deflection 
Xjj  is  approximately  equal  to  .iThat  it  iTOuld  be  in  the  lim ting  case, 

■  ■  ^II  “  s|  Po^o- 

According  to  Sq,  (21 )  ve  msh  to  ccxuparo  the  quantities 


'  ;'fPo)’ 


T)  ^  Po^o^i- 


h  ! 


P  \3' 


It  is  clear  that  in  general  these  ti/o  quantities  do  not  approximate  each  other, 

and  as  a  fiorther  check  it  is  easy  to  choose  special  values  of  the  constants 

v/hich  shcTw'  a  substantial  discrepancy  betireen  the  tvTO  quantities. 

I  PI  ^II 

As  a  numerical  example  suppose  that  M=l,  to '"I,  ~  -  and  ■::: - 100.  ‘ 

Po 


Then 


Po 
•II  12 » 

III 


X. 


X-rxT  =  hxi. 


and 


^=100 
Xi  12xi» 


Pq  =  I200x^; 


Vfe  msh  to  compare  the  quantities 


l6x|  and 


Po^i' 


or 


or 


16x1  and 


Po 


16x1  and  200xi. 

These  quantities  differ  by  a  factor  of  more  than  12.  For  thJ.s  numerical  case 
practically  all  of  the  action  takes  place  t.MIo  the  right-loand  side  of  the 
differential  equation  is  positive . 


rr 


OAtt 

U.t.  OAB 

COCNTIT 

lANOUACI 

pAoe 

RunruTioNs 

Auk  •  45 

Unclass. 

u.s. 

EnRllsh 

17 

tables,  graph 

TITLErAir  am)  Earth  Shock  -  A  Compilation  of  Informal  Reports  Submitted  in  Advance  of 
Final  Reports  -  Vol.  12  •  25  June  to  25  July  1945 
AUTHOR(S)  :  (Not  known! 

ORIG.  AGENCY  :  Princeton  Unlv.,  N.  J. 

PUBLISHED  BY  :OSRD,  NDRC,  Div.  2 ,  Washington,  D,  C. 


ATI-  60554 


•tVHIOM 

(None) 


OtIO.  ACOtCT  NO. 

(None! 


rwMSMIMO  NO. 


ABSTRACT; 


The  differential  equation  Mx  4^  F(x)-p(t)  is  considered  for  some  simple  cases  of  blast 
loading.  The  right-hand  side  is  assumed  linear,  and  F(x)  on  the  one  hand  is  taken  as 
constant  and  on  the  other  is  taken  as  linear  from  the  origin  to  the  constant  and  then  as 
remaining  constant  for  larger  values  of  x  .  It  is  shown  that  the  situation  In  the  two  cases 
differ  moderately.  An  approximation  formula  is  deTelq;)ed  by  udiich  certain  information 
in  the  latter  case  can  be  obtained  from  the  former. 


DISTRIBimON;  Copies  of  this  report  obtainable  from 


DIVISION:  Army  Materiel  (27) 

SECTION:  Army  Ordnance  (1) 

ATI  SHEET  NO.:  R-27-1-2 

SUBJECT  HEADINGS:  Explosions  -  9iock  and  pressxtre 
measurement  (34458);  Blast  effects  (16210);  Stock  waves, 

One -dimensional  (86309.25);  Pressure  waves  (74620) 

Control  Air  Docomontt  OMco 

1  Wrieitt<#nMonon  Air  Fereo  Bo«',  Dayton,  Ohio 

AIR  TECHNICAL  INDEX 

J 

